Pattern Formation in Phase-Separating Gels with Spontaneous Shear 



Nariya Uchida* 

Max- Planck- Institut fur Polymerforschung, Postfach 3148, D-55021 Mainz, Germany 
(Submitted August 7, 2001; published as Phys. Rev. Lett. 89, 025702 (2002)) 

We study pattern formation in gels undergoing simultaneous phase separation and orientational 
ordering. A 2D numerical simulation is performed using a minimal model of nonlinear elasticity 
with density-anisotropy coupling. For strong positive coupling, the collapsed phase elongates along 
the phase boundary and buckles, creating a folded structure with paired topological defects. For 
negative coupling, soft elasticity of the swollen phase causes a droplet morphology as in liquid-liquid 
phase separation. Their possible realizations in nematic liquid-crystalline gels are discussed. 

PACS numbers: 64.75.+g, 82.70.Gg, 83.10.Tv, 61.30.Jf 



Polymer gels are unique objects that undergo large 
deformation upon phase transition, which creates vari- 
ous patterns. Characteristic to pattern formation in gels 
is the interplay between rubber-elasticity and interfacial 
tension, as demonstrated by the surface undulations seen 
in swelling/shrinking processes Q. Elasticity also affects 
phase separation in bulk, which produces a foam-like 
morphology with the polymer-rich region connected [^|,^| . 
A similar morphology is observed in viscoelastic polymer 
solutions Q as well as in polymerization-induced phase 
separation [||. While these patterns are caused by vol- 
ume phase transition of the polymer network, there is 
a class of amorphous solids in which transition induces 
shear deformation. The resulting low-temperature phase 
is an intermediate between usual solid and liquid, in the 
sense that one of the linear shear moduli vanishes in that 
phase ||. We shall refer to it as "anisotropic glass" (AG) 
after Ref. ||. An example of AG is provided by ne- 
matic liquid-crystalline elastomers fj]||, which elongate 
along the director upon the isotropic-nematic transition. 
In their phase ordering process, the spontaneous shear 
deformation creates an anisotropic packing structure of 
orientationally correlated regions, which gives rise to an 
extremely soft nonlinear elastic response [p|jic|l. The 
anisotropy axis of deformation is locally parallel to the 
nematic director, which can be integrated out to give 
an effective free energy as a functional of strain. The 
aim of this Letter is to address pattern formation in 
phase-separating AGs, where volume change and spon- 
taneous shear cooperate plf . Motivated by the idea of 
cell-dynamical systems [[L2|, we adopt a minimal model 
of non-linear (visco-)elasticity that retains the essential 
features of (i) fixed-point configuration in the order- 
parameter space and (ii) vectorial volume-conserving dy- 
namics. 

We consider a two-dimensional system for simplicity. 
Let r and R be the positions of material element be- 
fore and after deformation, respectively. The field R(r) 
defines the deformation tensor dRi/drj and a symmet- 
ric strain tensor, GV, = (8Rk/dri)(dRk/drj). For a 
system with rotational invariance, free energy change 



due to homogeneous deformations can be expressed in 

1/2 

terms of the two invariants a = (det G) = A1A2 and 
t = trG/2 = {\\ + A|)/2, where Ai and A2 are the prin- 
cipal elongation ratios. Note that a is the volume expan- 
sion ratio while r — a = (\i — X 2 ) 2 /2 gives a measure 
of shape anisotropy. In dimensionlcss units, the homo- 
geneous part of our model free energy density is written 
as 

/ = (a d P f{d d Q f + ( Mo + r + ^r 2 . (1) 

The first term with the constants dp and ctq(> op) in- 
duces density phase separation, while the coefficients /io, 
fii, and vq control shear elasticity. When /io > and 
fix = vq = 0, the form of the shear free energy agrees with 
what prescribed by the classical affine-deformation the- 
ory of rubber-elasticity The coefficient /ii accounts 
for the coupling between density and shape anisotropy, 
and vq is introduced to ensure stability when /io or /ii is 
negative. 

Because of the nonlocal nature of elastic distortion, 
phase coexistence in a gel crucially depends on boundary 
conditions and geometry, and admits no simple and exact 
analysis like Maxwell construction fl4) . Here we restrict 
ourselves to the configuration of free energy minima in 
the a — r plane, which helps classification of phase behav- 
ior and drawing an approximate phase diagram. For this 
purpose it is useful to define the effective shear modulus, 

df 

V = q- = Mo + V\o + vqt. (2) 

Minima of / under the constraint a < r is controlled 
by the location of the line /i = with respect to the 
points P : {dp, dp) and Q : {dQ,dQ) on the d — r 
plane (Fig.l). There are four possibilities: (i) isotropic- 
isotropic (I-I) separation: P and Q locate the minima 
if both fi\p and /x|q are positive; (ii) AG-I separation: 
if fj,\p < and /i\q > 0, the free energy minimum cor- 
responding to the collapsed phase is on the line fj, = 
and off the isotropic line <r = r; (iii) I-AG separation 
(p\p > 0, /x|q < 0): the swollen phase is accompanied by 
spontaneous shape anisotropy; and (iv) AG- AG separa- 
tion (/i|p < 0, /i|<g < 0). 
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The total free energy is the sum of homogeneous and 
gradient contributions, and is assumed in the form 



F = dr 



/(a,r) + ^(W) 2 + y(VG) 2 



(3) 



While the second term penalizes density gradients and 
drives domain coarsening, the last term penalizes also 
gradients of the principal strain axis, and drives orien- 
tational ordering. To describe the latter, we define a 
"nematic" order parameter Qij = Gy — GkkSij/2 @- 

For the dynamics, we assume the simplest equation for 
gels with viscous solvents {ana] , 
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where T is the inverse of viscosity. Linearizing with re- 
spect to the displacement u = R — r, we have the initial 
growth of the longitudinal mode uii (q) = q~ 2 qq ■ u{q) as 
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U|| (q) = — r [V + W + (L + AM)q 2 } q 2 u\\ (q), (5) 
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while the transverse mode u±_ = u — uu grows as 
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These give the spinodal conditions for density separation 
and orientational ordering as V + W < and W < 0, 
respectively. 

To study the nonlinear regime, we implemented the 
model on a 256 x 256 square lattice with periodic bound- 
ary conditions. The kinetic equation (^) was time- 
integrated with the Eulcr scheme. The reference param- 
eters are a P = 0.6, a Q = 1.5, L = 0.2, M = 0.05, 
/Ho = —0.11, fix = i/o = 0.05, and T = 0.15, with the 
time increment and grid size set to unity. For this set 
of parameters, the initial state is unstable both to phase 
separation and orientational ordering (as V — —0.23 and 
W = —0.01). For any field quantity A(r), we shall de- 
note its spatial averages over the collapsed (a < 1) and 
swollen (a > 1) regions by {^4}< and {^4}>, respectively. 

Figure 2(a) shows evolution of the network density 
field l/a(r). At an early stage, there appears a foam- 
like pattern connected by strands of the collapsed phase. 
Then the strands start to undulate. At each vertex of 
the foam-like pattern, the strands fold onto themselves 
to form a lump embracing two spots of the swollen phase. 
Finally the pattern coarsens in a self-similar fashion, 
with small vertices and thin strands dissolving. Plot- 
ted in Fig. 2(b) is the distorted lattice mesh defined by 



R(me x + ne y )(m,n: integers), and the corresponding 
"Schlieren" texture Q 2 y - Note that each folded structure 
contains a pair of 1/2 defects. 

To see the origin of the morphology, next we study the 
crossover between AG-I and I-I separation. The domain 
boundary becomes less wiggly as fj,o is increased, and 
the effective shear modulus in the collapsed phase {/z}< 
turns positive at = —0.076; see Fig. 3. For larger val- 
ues of /io, the model reproduces the foam- like morphol- 
ogy known for I-I separation 0]. In the bottom of the 
figure, we plot versus no the principal elongation ratios 
in the two phases. In I-I separation, the swollen phase 
is almost isotropic (Ai ~ A2), while the collapsed phase 
has a sizable strain anisotropy. This anisotropy comes 
from the volume mismatch between the two phases. At 
a flat phase boundary, continuity of the network mesh 
requires elongation of the collapsed region parallel to the 
interface, which creates elastic tension; see Fig 4(a). 

In AG-I separation, on the other hand, the AG re- 
gion chooses its anisotropy axis along the interface to 
reduce the tension. The boundary stress is completely 
canceled when the ratio of the spontaneous elongation 
matches that of the boundary-induced anisotropy. This 
is roughly identified with the case {Ai}< = {Ai}>, which 
occurs at ^0 = —0.088. For larger anisotropy, or when 
{Ai}< > {Ai}>, the mismatch stress is canceled by tilt- 
ing the strain axis in the AG phase, as illustrated in 
Fig. 4(b). At boundaries between interfacial regions with 
different tilt directions, however, continuity of the net- 
work mesh requires a buckling of the interface. At the 
same time, torque is created at vertices where different 
strands meet, since the tilt is opposed by the clamping 
at the center; see Fig 4(c). This torque causes a spon- 
taneous winding of the structure, which is further cut 
at two points near its center to reduce elastic free en- 
ergy, creating two defects (Fig.4(d)). This leads to the 
patterns as seen in Fig 2.(b). 

Now we consider the case of I- AG separation by as- 
suming a negative density-anisotropy coupling. Shown 
in Fig. 5(a) is the coarsening process for the parameters 
fio = 0.05 and fi± — —0.11. Initially a foam-like mor- 
phology like in I-I separation appears. As anisotropy 
in the AG phase develops, the cell walls break and the 
swollen domains merge. Finally the isotropic phase form 
droplets, as is the case for liquid-liquid phase separation 
with the same volume fraction. In Fig.5(b), we illus- 
trate the deformation induced around an droplet of the 
isotropic phase. To reduce the mismatch at phase bound- 
ary, the principal elongation axis in the AG phase tends 
to be perpendicular to the interface. Therefore, the de- 
formation in the AG phase is necessarily inhomogeneous. 

The formation of the droplet morphology is understood 
from the fact that an AG in inhomogeneous deformation 
states is liquid-like soft |T^,|l6|. Internal elastic stresses 
are greatly reduced due to a screening of strain-mediated 
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long-range interaction between orientationally correlated 
regions. In the present case, it means that stresses cre- 
ated at a phase boundary are not transmitted far into 
the AG regions. Thus, the foam-like pattern cannot be 
mechanically supported and the morphology is controlled 
by interfacial tension, leading to the liquid-like pattern. 

To summarize, novel patterns should arise in phase- 
separating gels with density-anisotropy coupling. To 
realize the wiggly self-folded morphology, we need the 
anisotropy ratio A1/A2 to be 3 ~ 4 in the AG phase. 
It is a realistic value in nematic gels made up of main- 
chain liquid-crystalline polymers jL7|. Also, negative cou- 
pling and collapsed-isotropic-swollen-nematic separation 
can be achieved by using nematic solvents fig}] . The 2D 
patterns should be observable in a film geometry as used 
in fabrication of polymer-dispersed liquid crystals, which 
basically are isotropic gels with nematic solvents and do 
show quasi-2D patterns fl. It should be interesting to 
extend the study to a 3D system. There the swollen 
isotropic phase creates an isotropic tension along a flat 
phase boundary, while uniaxial deformation of the AG 
phase can cancel it only in one direction; the stress in 
the other direction would cause a spontaneous curvature 
of the boundary. On the other hand, the liquid-like mor- 
phology would be qualitatively unchanged, since the soft 
elasticity of the AG phase is due to the presence of bulk 
Goldstone modes. 
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FIG. 1. Schematic configurations of free energy minima 
(marked by •). (a) I-I separation, (b) AG-I separation, (c) 
I-AG separation, (d) AG-I separation with the initial state 
stable against orientational ordering (see Eqs.(7) and (8)). 
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FIG. 2. (a) Evolution of the network density field. Col- 
lapsed (AG) regions are plotted in black. Shown is a 192 x 192 
portion of the distorted lattice, at t = (1, 4, 16, and 64) x 10 4 
from the top left to bottom right, (b) Distorted lattice mesh 
(red) and "Schlieren" texture Q 2 xy (grayscale). 
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FIG. 3. Crossover from AG-I to I-I separation. Top: do- 
main morphologies at /io = —0.09 and fio = —0.05. Middle: 
the effective shear modulus averaged over collapsed regions, 
{/i}<. Bottom: principal elongation ratios averaged over 
collapsed and swollen regions. All the data are obtained at 
t = 32 x 10 4 . 
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FIG. 4. Deformation at interfaces, (a) If the preferred 

elongation ratio Ai of the collapsed phase is smaller than H _ l°J~l" ~ r*i~*". t I ~ 

that of the swollen phase, mismatch of the network mesh 

creates tension along the interface, (b) In the opposite case, 

the mismatch stress is canceled by tilting the elongation axis ' ' | ' 

in the AG phase, (c) Torque is created at each vertex, where 

tilted strands are clamped, (d) The cross-shaped structure 

in (c) winds up itself and is cut at the points marked by o. 11111 

FIG. 5. (a) Density snapshots for I-AG separation, at 
t = (1, 4, 16, and 64) x 10 4 from the top left to bottom right, 
(b) Boundary matching induces inhomogeneous anisotropy 
in the AG phase around an isotropic droplet. 
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